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I. INTRODUCTION 



In hydrography, the character of the bottom should be determined for nautical 
charting, particularly in harbors, anchorage areas or ether areas where the bottom 
characteristic is a significant factor to be considered for navigation safety. 

Conventional samipling methods, using clamshell snappers (Umbach, 1976). are used to 
define the characteristics of the bottom surface layer. It is tedious and time 
consuming, especially in oIT-shore geologic investigations. If some kind of theory or 
model existed to remotely classify the sediments, the task could be carried out easier 
and faster. .Moreover, this ability would permit nautical charts to contain much more 
information. 

In the ocean bottom, the magnitude and nature of the sound scattering are 
functions of both the particle size of the bottom and surface (bottom) relief (.Mckinney 
and Anderson. 296-j. The sound eclio from sediments is subject to tliree types of 
reflecting or scattering (Clark, Proni. Seem, Tsai. 1985) : 

1) coherent reflection from the water-sediment interface; the signal amplitude 
will be approximately that given by spherical spreading of the wave front 
from the image source modified by the plane wave reflection cocflicient 
(Clay, 1966), 

2) incoherent or statistically variable scattering from the surface irregularities 
of the sediment, and 



3) incoherent scattering from within the volume of the sediment caused by the 
acoustically irregular matrix of the sediment. 

Items 1 and 2 arise from the surface of the sediments and item 3 from the sound 
pulse penetrating into the sediment volume and then being scattered back to the 
transducer by the granular matrix of the sediment. Since different sediments have 
different acoustic matrices, they have different volume reverberation characteristics and 
there should exist a mathematical relationship between sediment properties and the 
volume reverberation. 

The major goal of this experiment was to establish a fundamental laboratory 
procedure for the remote sensing of ocean sediment volume reverberation. This was to 
be accomplished by accurately measuring the properties of a pulse reflected from the 
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sediment. In application, since the I'dlected sound pulse is received by the same 
transducer, i.e., "monostatic sonar", our discussion is restricted to "back-reflection" and 
"backscattering". 
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II. BACKGROUND 



Bradshaw (1981) used Monterey #31 fine sand, the same material used in this 
experiment, to study the sound propagation into a fast bottom medium. In hi.s study, 
air bubbles were removed from within the sediment based on heating the water sand 
system. Bleach was added at a ratio of one-half gallon bleach to 70 gallons of water 
plus sand in order to control biologic growth. The physical properties, such as density, 
porosity, speed of sound, rellecton coefficient etc., were measured. Kosnik (1984) used 
a high pressure water jet to remove the air bubbles from the water, sand mixture. 
However, no volume scattering measurements were made by either investigator. 

Morse and Ingard (1968) give a theoretical derivation on incoherent scattering. 
Under the assumptions that : 1) the scattering objects are quite small compared with 
wave length X. 2) the dimensions of scattering region are much larger than X, 3 ) the 
scatterers are all spheres, and 4) the scatterers are populated sparsely and no mutiple 
scattering existed, the incoherently scattered intensity Ij^. per unit incident intensity Ij. 
at angle 8 to the incident wave, for frequency co, 2rt, is 



where N is the total number of scatterers within unit volume, K is the wave number for 
the mJxture, a is the radius of scatterers. Furthermore, for monostatic sonar, the 
volume scattering coefficient for unit volume can be reduced as 




( 2 . 1 ) 



(J_K^ a" sin"( 1, 2)0) exp(-2K“ a" sin"(l '2)8) 
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( 2 . 2 ) 



where 



K 



n 



K 



0 



K 



0 
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3 Pn - 3 Po 

Yp = 

2 Pn - Po 

and is the compressibility of the scatterer. Kq is the compressibility of the fliiid. 
is the densitv of the scatterer and is the densitv of fluid. 

For a small (ka < < 1) rigid sphere, the backscattering cross section is given 
{Urick, 19S3) as 

a = 2.S (It a2){K a)^ 

Note that, as a single sphere, it must be a single scattering. The scattering coenicient 
for unit volume can be obtained as 

s^. = N o (2.3) 

Both theories are not very appropriate for the real ocean sediment volume 
reverberation. However, these are the only two theories available. The backscattering 
of sound from the ocean bottom has' been studied for many years, but unfortunately no 
papers were found to discuss the volume reverberation of ocean sediments. 
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III. DESCRIPTION OF VOLUME REVERBERATION 



As sound propagates in a fluid, part of the sound is intercepted and reradiated by 
inhomogeneites in the fluid. The reradiation of sound is called scattering. The total of 
the scattering contributions from all the scattercrs is called reverberation lUrick. 1983). 

Volume reverberation is defined as sound scattered back to the transducer by 
scattering centers in the volume of the sea (National Defense Research Commitee, 
1969). Volume reverberation may be considered as a blending of large number of 
echoes. It has been studied mostly in the Deep Scattering layer . a region of high 
biological activity. Reverberation level, in decibels, is used to measure volume 
reverberation. Consider a directional projector in an ocean containing a large number 
of volume scatterers trigtite 3.1). Let the axial intensity at unit distance be 1^. At a 
distance r where there is a volume element dV of volume scatterers, the incident 
intensity is 

Jr = lo 

The scattered intensity 1 m from the volume dV is 
dls(l)= li-S-dV 

where s,. is the volume scatterin 2 coefficient. The scattered intensity back at the 
receiver is 

dl^ (r) = Ij s,. dV;r2 

= (Ior^)s,.dV 

and the scattered intensity from all scatterers contributing to the scattered intensity at 
time t is 

I, (r) = s,, V 

where V is the total volume contributing for the geometry of Figure 3.1 calculated by 
(L’rick, 1983) 

V = F2 r- ct;2 
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Figure 3.1 Portion of ocean scattering. 



IS 



where is the solid angle of the sound beam. T is the pulse length and c is the sound 
speed, so 



^s (^0 ^ ^ 

The reverberation level is 

RL = 10 LOG Is(r) lj.(,p 

= SL - TL + S^. + 10 LOG (H ct 2) 



where TL is the transmission loss for spherical spreading 20LOG(r) and is the 
scattering strength for unit volume 10LOG(s^.). In Deep Scattering Layer studies, this 
reverberation level equation is sulTicient, because the received pulse is regarded as 
totally coming from volume reverberation, i.e. negligible reflection and surface 
scattering. Furthermore, the derivation of V is based on a square acoustic pulse. But 
for a real transducer, the square pulse is transformed into a pulse with an exponential 
rise, a flattened top and an exponential decay. The sound pulse reflected from the 
sediment and received by the transducer combines coherent reflected, incoherent 
scattering from the surface and volume reverberation. 

Let a transducer located in the water send out a sound pulse to a smooth 
water air interface. The received echo waveform is purely rellection. If the water 
interface is rough, then the received echo waveform includes both reflection and surface 
scattering. Since sound energy is scattered by the rough surface, the received waveform 
amplitude is reduced. Consider first-order scattering and suppose that the surface 
roughness is the same for every position. The amplitude of the combined sound wave 
(reflected and scattered) is reduced proportional to the purely reflected wave amplitude. 
In Figure 3.2(a). the solid line represents an ensemble of reflected waveform, and the 
dashed line represents an ensemble of reflected and scattered waveform. 

The contributions to a pulse reflected from a sediment surface is shown on 
Figure 3.2(b). The dashed line is the contribution by reflection and surface scattering; 
again, only first order scattering is considered. The dotted line includes volume 
reverberation. The level of the signal is raised above that of a dashed line because of 
energy returned from within the sediment. In particular, the tail decays more slowly 
than for reflection from the water'air interface. This is the significant feature of 
sediment volume reverberation that we are after in this paper. 
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Figure 3.2 '1 heoretical echo waveform for water and sediment. 



20 



Two alternative assumptions have to be made for reflection from the sediment 
surface : either !) the surface scattering is negligible compared with volume 
reverberation, or 2) if it is not negligible, it must have the same amplitude factor as 
discussed in Figure 3.2(b). 

According to these assumptions, if a sediment echo is subtracted from the water 
echo it will give a measure of the volume reverberation in the sediment (Figure 3.3 i. 

The decaying tail of the echo reP.ected from the smooth water, air interface is an 
exponential decay. 

V„. („ - 

where (t) is the received voltage at time t. 6 is the decay constant. is a constant 
voltage, and t is the time starting from the beginning of decay. If we assume the 
sediment echo also decays exponentially, then we will obtain a similar function 

where (t) is the received voltage from sediment at time t and J] is a decay constant. 
If and Vj are both normalized to the same initial value K. the dilTerence becomes 

\\. (t) = K (e-'l^ - e’^^) (3.1) 



If both Vjt and 5t are much less than unity. 

\\. (t) = K ( n - 5 )t t < < 1 J) and 1 6 
This initial decay is linear and depends on the dilTerence 6-i]. 
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ECHO FROM WATER SURFACE 



I 





Figure 3.3 Real echo waveform displayed on CRT. 
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IV. EXPERIMENT DESIGN 



A. MATERIAL AND TANK SELECTED 

Fresh (tap) water, t?30 Monterey fine sand and Monterey Aquarium ^2 for a 
rough gravel sediment were the media used in the experiment. The grain size of ==30 
fine sand varies from 0.15 mm to 0.7 mm (Bradshaw, 1981) with an average size 0.3 
mm. The Aquarium ^2 has a mean size 5.3 mm (Diaz. 1986). The sand water inixturc 
was repeatedly stirred to remove air bubble trapped while transferring the sand to 
different tanks. 

Initial experiments were conducted in the NTS anechoic tank to measure the 
reflection from the water/air surface, with the transducer located at dilTerent depth. 

The second tank, a wooden tank measuring 80 cm x 80 cm x 60 cm, was constructed to 
hold about 30 cm of aggregate. This tank was lowered and suspended 85 cm below the 
water, air surface of the anechoic tank. A third tank, a steel glass tank, measuring 70 
cm X 70 cm x 60 cm was filled with 24 cm of fine sand or aggregate with 36 cm of 
water above. 

B. ELECTRONIC EQUIPMENT 

A schematic drawing of the equipment is shown in Figure 4.1. E.xcept a high 
pass filter, all components were off-the-shelf General Radio model 1310 oscillator 
with a frequency of 182.0 kHz was fed simultaneously into a Hewlett-Packard 5233L 
frequency counter, and a General Radio Type 396-A tone burst generator was used to 
generate a 16 cycle pulse. Signal output from Hewlett-Packard HP 467-A amplifier 
was fed through a Datasonic Transmit Receive (T;R) switch and a small 21 to 28 kHz 
high pass filter, then into the F-41 transducer. The F-41 transducer was discussed in 
detail by both Borchardt (1985) and Diaz (1986). Signals returned from the sediment 
were amplified 20 dB by a Hewlett-Packard 465-A pre-amplifier, then passed through a 
Spencer- Kennedy Laboratories, Inc. model 302 variable electronic filter (set at 135 kHz 
high pass) to eliminate low frequency noise before being passed to the digital and 
analog oscilloscopes. 

The waveform displayed on 309 1 Nicolet Digital Oscilloscope was sampled and 
sent as a stream of characters to the computer. To reduce manual operations on the 
Nicolet 3091 during the experiment, a IIP-3421;\ Data Acquisition and Control Unit 
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Figure 4. 1 Schematic of equipment. 
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was adopted to receive commands from a HP-S6 computer and then automatically 
operate the Nicolet 3U91. However, it was noticed that if the Nicolet 3091 was not 
"warmed up" enough, the DC-offset jumped up and down, and continual manual 
adjustment was needed. .After 2 or 3 hours of running, it became stable and no further 
adjustment was required. 



C. VERIFICATION OF MEASUREMENTS 

Xo matter how deep the transducer is located below the water air surface, as 
long as the transducer is level, the exponential decay of the tail should be the same. In 
Diaz's thesis, the measurements were inconclusive because the transducer was fixed at 



one position throughout the experiment. This was not adequate to prove the system 
functioned properly. To prove that the previous experiments were run properly and 
results were obtained correctly, a verification was attempted by measuring the 
rellection from the water air interface with the same procedure, but with the transducer 
at various depths. Data set 20 (Figure 4.2) for tiie transducer 70 cm below the surface, 
gives a slope of-3.~5 x 10'“ Xp ;is. Data set 2! (Figure 4.3) was obtained by raising 
the transducer 15 cm and the slope is -3.77 x 10‘“ Xp }is. Both results were close to 
the previous results obtained by Diaz (1986). After these two data sets were collected, 
a new transducer mounting was built and adopted to improve the experimental 
precision and operation ability. The transducer was leveled by a counter- weight on the 
new carriage. To ensure the system and operation worked properly, more locations 
were tried to obtain the decay constant of the water, 'air surface. 

Data sets 23 and 24 were obtained at the same transducer depth, but seperated 
40 cm horizontally. Data sets 24 and 25 were obtained at the same horizontal position 
but with 50 cm difference in depth. Data set 26 was acquired by changing 10 cm 'in 
both the depth and position from data set 25. Data sets 26 to 28 resulted from 
changing the depth from 70 cm to 30 cm below the water air interface. The results of 
these experiments are listed in Table 1 and show consistencies between each other. 
Waveform ensembles of the tail for each data set are plotted on Figure 4.4 to Figure 
4.9. 



The average slope -3.77 ± 0.077 x 10’“ Xp fis is smaller than -3.84 ± 0. 14 x 
10’“ Xp. Jis which was measured in Diaz's thesis. Statistically, these two slopes are the 
same. The smaller standard deviation might suggest that the new carriage has achieved 
a better precision. 
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Figure 4.2 Decay envelope for icdeclion (Vom water/air interface, 
transducer 70 cm below the suiT'ace (data set 20). 
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Figure 4.3 Dccav envelope for lellection IVom water/air inleiTace. 
transducer 55 cm below tlie sui l'ace (data set 21 ). 
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2-OOr CORRELATION C0EFF- -.9795 



TABLE 1 

WATER AIR REFLECTION 



DATA 


SLOPE 


CORRELATION 


SET 


(xlO-2) 


COEFFICIENT 


DATA20 


-3. 75 


-0. 9984 


DATA2 1 


-3. 77 


-0. 9795 


DATA22 


-3. 86 


-0. 9979 


DATA23 


~3. 62 


-0. 9724 


DATA24 


-3. 68 


-0. 9895 


DATA25 


-3. 86 


-0. 9956 


DATA2 6 


-3. 80 


-0. 9679 


DATA2 7 


-3. 84 


-0. 9888 


DATA28 


-3. 79 


-0. 9810 


MEAN SLOPE = 


-3. 77 X 10“2 Np/Hs 





a = 0. 07 X 10"^ Np/^s 



After dealing with the water air surface, the transducer was turned over and 
pointed downward toward the sediment which was contained within the 80 cm x SO cm 
X 60 cm wooden tray. Difierent transducer positions were choosen. Data sets 31 to 33 
were at the same depth, but horizontally 3 cm away from each other. The results show 
the waveform were totally inconsistant. Data sets 33 to 35 were at the same horizontal 
position but dilTerent depths; results were also inconsistant. Above results are 
presented on Table 2 and plotted on Figure 4.10 through Figure 4.14. So far, 
experiments reveal the lack of isotropic homogenicty of the volume contributes 
reverberation. 

D. MEASUREMENT CONSIDERATION 

During the course of these experiments some other effects were considered and 
studied to ensure that the experiment was working properly. 
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Figure 4.4 Dccav envelope for rellcction lioni water/air interface, 
transducer was 30 cm below surface (data set 23). 



29 



CORRELATION C0EFF= -.9724 



TIME (MICRO-SEC) 



LN((<VOLTAGE'2>) -1/2) 




Figure 4.5 Dccav envelope for leflcction from water/air intciTacc, 
trnsduccr was 3U cm below surface (data set 24). 
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CORRELATION C0EFF= -.9895 




Figure 4.6 Decay envelope for rcllection from watcr/air interface, 
transducer was SO cm below suilace (data set 25). 
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I''igure 4.7 Decay envelope for rellection I'rom water/air interface, 
transducer was 7() cm below surface (data set 26). 
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2-OOr CORRELATION C0EFF= -.9679 
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Figure 4.8 Decav envelope lor rdlcction I'roni \vater;air intciFacc, 
transducer was 50 cm below suri'ace (data set 27). 
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2. OOr CORRELATION COEFF= -.9888 




Figure 4.9 Dccav envelope for rellcction from watcr/air interface, 
transducer was 3U cm below surface (data set 28). 
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TABLE 2 

SEDIMENT WATER REFLECTION AND SCATTERING 



DATA 


SLOPE 


CORRELATION 


SET 


(xl0“2) 


COEFFICIENT 


DATA3 1 


-2. 20 


-0. 9138 


DATA32 


-2. 98 


-0. 8908 


DAT A3 3 


1 — 1 
1 


-0. 7838 


DAT A3 4 


1 

o 

0^ 


-0. 8174 


DATA35 


-2.90 


-0. 9844 


MEAN SLOPE = -2.62 : 


t 10“2 Mp/ps 





(7 = 0. 86 X 10"^ Np/us 

1. Effects of Sampling Rate 

The sound signal being used was fixed at 1S2 kllz (5.5 fis period). The 
Nicolet 3091 digital oscilloscope has as its fastest sampling rate 1 ps which is 
inadequate to give a precise reproduction of the numeric waveform. However, 
observations have shown "jitter" in the triggering of the scope, so for consecutive 
pulses the scope would sample dilTerent parts of the waveform. To obtain an accurate 
waveform, a large number of waveform were sampled and averaged. If sufiicient 
"jitter" is present, the voltage measured in a given bin should van; between the 
maximum and minimum voltage in the waveform in the vicinity of the corresponding 
time. Since the sine wave varies between minus and plus, the voltage values were • 
squared before averaging in order to give an accurate representation of the waveform 
envelope. 

To determine the maximum number of samples required for each data set 
(average waveform), the water'air interface was measured with different number of 
samples. The results showed that 50 samples for each data set were enough to give a 
true waveform. 

2. Effects of Bottom Reflection Interference 

For a 16 cycle signal with a frequency 182 kllz, the pulse length is 88 )ls. 

Since the sound pulse emitted had a decay time of appro.ximately 100 fts, the echo from 
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Figure 4.10 Dccav envelope for rellection from aesiegate, 'water 
interlace, transducer was 30 cm above surface (data set 31). 
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Figure 4.1 1 Decav envelope for rellection from aegregatcAvater 
interface, transducer was 30 cm above surface (data set 32). 
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2-OOr CORRELATION COEFF- -.8908 
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Figure 4.12 Decay envelope for rellection from aegregate, 'water 
Intel (ace, transducer was 30 cm above surlace (dfita'set 33). 
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Figure 4.13 Dccav envelope for rellection from aggregate/water 
'interlace, transducer was .50 cni above surface (data set 34). 
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Figure ‘1.14 Dccav envclojic for rcllcclion IVoni aggregate, 'water 
interface, transducer was 70 cm above surface (uata'set 35). 
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2. OOr CORRELATION C0EFF= -.9844 



the smooth water air interface should also have the same decay time. Since the 
purpose of this project was to compare the decay portion of echos rellected from water 
and sediments, the waveform was sampled for a 100 }ts interval from the begining of 
decay. The sound speed in the aggregate was measured to be 1555 m s. For the fine 
sand. Bradshaw (19S1) obtained the speed of sound to be 1610 m s. Combining t'nese 
sound speeds and the time interval ISS us. the depth of sediment that contrubutcs to 
the main echo is 14.6 cm for aggregate and 15.1 cm for fine sand. Interference caused 
by reflection from the bottom of the tank will not alTect the experiment as long as the 
sand is thicker than 15.1 cm and the aggregate is thicker than 14.6 cm. Throughout 
our experiments, we kept both sediments more than 24 cm deep in the glass tank to 
avoid interference. 

3. Effects of Near Field Interference 

In the near field of a transmitter, the axial pressure exhibits strong interference 
eflects (Kinsler. Frey. Coopens. Sanders. 19S2). Thus, the sound field is irregular and 



does not fall off smoothly with distance as it does in the far field. To avoid 
measurements in the near field, an experiment was run to find the transition distance 
which separates the near field and far field. I'his distance is 20 cm from the transducer 
center. No problems should occur as long as the target is more than 20 cm from the 
transducer. In the glass tank, we placed the transducer 30 cm above sediment. The 
water air interface measurement also gave a check that 30 cm distance is satisfied. 

Since a pulse moving in water at 14S0 m s, takes 405 ps to travel the 30 cm to the 
sediment interface and back, the 30 cm distance is long enough so that the transducer, 
which is on about ISS ps, stops ringing before the first echo from the sediment arrives. 

4. Effects of Water/ Air Interface Interference 

Ideally, the transducer should be far enough away from all unwanted surfaces 
to avoid interference. As the transducer faces down and projects a major part of its 
energy to the sediment surface, some energy is transmitted to the water surface and 
reflected back toward the sediment. If the transducer is close to the water, air surface, 
the reflected energy will interfere with the transmitted energy. Restricted by the size of 
the glass tank, we had to locate the transducer close to the water surface. To test how 
important this interference could be, the water surface was agitated by fingers. The 
waveform displayed on both analog and digital oscilloscope was unchanged during the 
ruffling. It was concluded that this source of interference is negligible for this 
experiment. 
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5. Effects of Side Lobe 

The half beamwidth (to the first null) for the F-41 transducer was measured to 
be 10.0’. For a 30 cm distance from sediment to the transducer, the sediment surface 
ensonified by the major lobe is a circle area with radius 5.3 cm. If the distance is SO 
cm. the radius will be 14.1 cm. For both the glass tank and wood rank, if the 
transducer is located above the center of the tank, the major lobe will not produce 
interference caused by the sides of the tank. 

The glass tank was too small to avoid the interference, if any, caused by the 
side lobe. To test how much the interference was. sediments outside the center 30 cm \ 
30 cm square area were disturbed and left in a ven.' irregular surface. Waveforms 
returned from sediment were measured before and after the disturbance. Data set 51 
(Figure 4.15) shows the waveform measured before the sediment was disturbed and 
data set 52 (Figure 4.16) shows the waveform after disturbance. By comparing the two 
waveforms, it is found that the F-41 trandsucer gave negligible side lobe interference in 
this experiment. 

6. Effects of Sediment Inhomogeneity 

Before the test findings on side lobe effects, we were bothered by the drastic 
changes on the waveform whenever the transducer was moved parallel to the smoothed 
aggregate surface. It was assumed that the waveform fluctuation was caused by the 
inhomogeneitics of the aggregate. To prove this, we placed the transducer at a fixed 
position in the glass tank and repeatedly ensonified the aggregate. After each data set 
was acquired, the aggregate was disturbed and smoothed again. Ten sets of waveform 
data were obtained and each of them differed from the others. 

It was thought the fine sand would be more homogeneous than aggregate and 
it might produce more consistant waveforms as the transducer was moved horizontally. 
Measurements made with the fine sand, however, showed that the waveform still 
fluctuated. After each data set was obtained from the smoothed fine sand, the 
transducer was moved horizontally to different position. .All the transducer positions 
were at least 15 cm from the wall to avoid interference caused by reflection from the 
sides of tank. Also, the positions were at least a length of tranducer diameter from 
each other to avoid a waveform replication. .A metal plate was set on top of the sand 
and the transducer was moved horizontally. As expected, the waveform on the 
oscilloscope was unchanged but totally different from sand. 
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Figure 4.15 Decav envelope for redection from aggrcgale/watcr 
iiuerlace in the glass tank (data det 51)7 
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Figure 4.16 Decav envelope for redection from aegregate, 'water 
interlace in the gla“;s tank (data set 52). 
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CORRELATION COEFF» .0400 



E. PROGRAMMING 



An HP-S6 desktop computer was adopted to collect and analyze the waveform 
data. The computer program was written to be interactive. This enables 1) the 
collection of data sets with din'erent delay time intervals and 2) the collecting and 
processing of numerous data sets. The interactive programs allowed essential 
infcrntation to be entered via the key board, eliminating the need to modify the 
program whenever data sets were changed. The interactive prompts were useful for 
people unfamiliar with Basic prograntming. 

1. Program "THESIS3" 

This program was designed as the initial step in the data collection processing. 
Before the selected waveform data displayed on the 3091 Xicolet Digital Oscilloscope 
were output to the computer, it was switched into "freeze" or "store" mode from "live" 
mode. Through the interface RS-232. the computer received the stored waveform data 
and then processed them. After processing the waveform data and before receiving the 
next sampling data, the Nicolet 3091 had to be switched into "live" mode again. An 
HP-3421A Data acquisition and Control L'nit was interfaced between the Nicolet 3091 
and HP-S6. The HP-3421.A reduced the amount of manual operations in the data 
collection process. A five volts source were input to the MP-342IA. Commands form 
the HP-S6 and 5 volts signal were used to control, by closing and opening the channel 
in the HP-3421A, the "mode switching" of the oscilloscope. To avoid errors, the 
oscilloscope should be warmed up enough while entering the inital information. Once 
this initial information was entered and the oscilloscope was warmed up enough: the 
whole system runs automatically. 

Waveform data output from the oscilloscope were a stream of characters. To 
obtain the voltage and time value, all these characters needed to be normalized and 
converted by the computer. A program, listed in the Nicolet 3091 operation manual, 
was used for data value normalization. A data set, the waveform of the decay portion, 
was obtained by squaring and averaging 50 samples of waveform data for each 
transducer position. The running average formula used in the program was 

< V2>^ = ((n-l)< V2 

where < V“ > is n samples averaged squared-voltage, < is n-1 samples 

averaged squared-voltage and is squared voltage for the nth sample. Note that 
the average window for the waveform is 100 fis from the beginning of the decay. The 
above formula was processed for each corresponding time bin. 
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During each data set processing, some error messages, such as STRING OVF 
(which stands for the string variable overflow), were displayed on the HP-S6 screen. 

This error messaged could be fatal to the data set if the data has not been stored on 
the disk. To avoid this problem, the program was designed so that after ever; 10 
samples processed the updated data were stored oit the disk over the previous updated 
data. The programi eliminated the problem of the disk being full before the final data 
set was stored on the disk. This prevention could increase the processing time, but 
compared to the time wasted from the above problems, it is worthwhile. The final 
result for this program was 3 data files : one stored the squared voltage, one for tim.e 
data and one for the number count of the samples. 

2. Program "MAXIMUM" 

In order to obtained the envelope of the waveform, this program was designed 
to pick up the local maximum of the voltage data. Tire voltages, stored in the data file 
from 'TI1ESIS3", were read and compared with the neighboring voltage data to decide 
the local maximum. The deternuned local maximums and corresponding time bins 
were stored in seperated files and printed out. The data files were difierentiated by 
number according to the order they were acquired. 

3. Program "PLOT3" 

After the local maximums were determined, the decay portion of the waveform 
can be plotted by program "PLOT3". The experiments, for this project, took place 
mainly in the glass tank. A decay envelope measured in the glass tank has an interval 
between 520 and 619 fis starting from transmitting time. Some exceptions existed and 
gave difTerent time frames. 

To determine the decay constant, the voltage data were converted 
logarithmically and plotted vs. the time 

V = V, e-«^ 

LN ( V ) = -at + LN ( ) 

If the cursor on the oscilloscope was not centered, the resulting envelope showed a 
"zig-zag" instead of a straight line for a perfect exponential decay. 

The least-square method was adopted to calculate the best fit of straight line 

y = ax +b 
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where a, which is similar to the decay constant, is the slope and b is the intercept on 
the y*axis. According to Beers (1953), 



a = - 



- 1 I >n 



b =- 



I Xn' I Yn - I (^n>n) 

- (I X, )2 



where k is the total number of pairs of (x.y) values, n is the number between 1 and k, 
y Xj^ is the sum ot the total x value, T is the sum of the total y value. V x^^" is the 
sum of the squared x value, y is the sum of the squared y value and y (■'^nYn'* 
the sum of the multiples of x^^ and y^^. 

Note that y x^^” and y (x^^)“ are not the same, y Xj^~ implies that each 
value of X is squared and then a sum is made of these squares, while y (x^^)^ implies 
that all of the values of x added together and then this sum is squared. There is an 
analogous distinction between y and y x^^y y^^. In our case, the x-coordinate 

is the time value and the y-coordinate is the natural-log of the voltage. To check how 
close the observed and least-square estimated y-valtie, compare, a correlation 
coefilcient was calculated (Chatfield, 1984) 



( Xn - X X yj, - y ) 

r 

(I ( ^n ■ ^ ( >n - >’ ^ 

where r is the correlation coeHlcient, x is the mean value of x and y is the mean value 
of y. Slope and correlation coelTicient were both printed with the graph. 

4. Program "AVERAGE" 

The received echo waveform varies from position to position because of the 
inhomogeneous sediment. To find the decay constant, we have to average out the 
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inhomogeneities to obtain a mean value of the waveform. This program performs a 
running average on the original data sets, producing the final averaged waveform of the 
sediment. 

The above four programs are given in Appendix A. 
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V. RESULTS 



The data sets for reflection from aggregate vrater interface are data 42 to dl. 
Figure C.l to Figure C.9, and Figure 4.15. Table 3 shows the slope obtained for each 
data set varied between data sets. To get a mean value for the slope, tlie running 
average method was used to process the data sets. The running average meilioJ also 
shows how fast the average processing became stable. Data sets D.A.T2 to DATIO, 
Figure C.IO to Figure C.IS, are the results .f'or each average process. The slopes are 
listed on Table 4. It is evident that the slope becomes stable after 8 sets of data are 
averaged. The final averged slope is 2.13 10'“ Np jis and is quite different from the 

value reflected from the water air interface. The difference between the two slopes 
resulted from, the effects of volume reverberation in the aggregate, and the scattering is 
sisnificant. 



TABLE 3 

ORIGINAL AGGREG.ATE REFLECTION 



DATA 


SLOPE 


CORRELATION 


SET 


(xlO"^) 


COEFFICIENT 


DATA42 


-0. 94 


-0. 5273 


DATA43 


-1. 35 


-0. 5874 


DATA44 


-2. 68 


-0. 7848 


DATA45 


-2. 64 


-0. 8509 


DATA4& 


-0. 12 


-0. 1162 


DATA47 


-1. 81 


-0. 3971 


DATA48 


-0. 95 


-0. 3990 


DATA49 


-3. 77 


-0. 8654 


DATA 50 


-2. 47 


-0. 9756 


DATA51 


-0. 16 


-0. 0584 



The data sets for reflection from fine-sand water interface are data 55 to 64. 
Figure D.l to Figure D.IO. The slope of each data set is listed on Table 5. Evidently. 
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TABLE 4 

AVEFLAGED AGGREGATE REFLECTION 



DATA 


SLOPE 


CORRELATI' 


SET 


{ xl0"2 ) 


CCEEEICi: 


DAT2 


-1. 19 


-0. 6051 


DATS 


-1. 97 


-0. 9720 


DAT4 


-2. 32 


-0. 9175 


DATS 


-2. 18 


-0. 9343 


DAT 6 


-2. 10 


-0. 9914 


DAT 7 


-1. 98 


-0. 9952 


DATS 


-2. 16 


-0. 9932 


DAT 9 


-2. 20 


-0. 9943 


DAT 10 


-2. 13 


-0. 9958 



the waveform changes very much for dilTerent data sets. The same running average 
method was used to remove the inhomogenities. Data sets DAT 12 to D.\T20 were 
produced from the average process. They are listed on Table 6 and plotted from 
Figure D.l 1 to Figure D.19. The slope became stable after 6 data sets were averaged. 
This might suggest that the fine sand has a greater homogeneity. The averaged slope 
for the fine sand is 2.37 x 10’^ Np fls and is greater than the aggregate. However, the 
slope difi'crence between them is not as distinguished when contpared to the slope 
dilLerence between the sediment and the water; air surface. 

According to Equation 3.1, the decay portion volume scattering, represented in 
voltage, for the aggregate is 

V, = K(e-21300t . ,-37700t ) 

and for the line sand is 

V = K { - e'37700t ^ 

The above relations only existed without surface scattering or the surface scattering 
was negligible. 
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TABLE 5 

ORIGINAL FINE SAND REFLECTION 



DATA 


SLOPE 


CORRSLATD 


SET 


( 


COEFFICi: 


DATA 5 5 


-1. 51 


-0, 8178 


DATA 5 6 


-3. 03 


>0. 9581 


DATA 5 7 


-2, 50 


-0. 8112 


DATA58 


-0. 75 


-0. 7290 


DATA59 


-1. 00 


-0. 6776 


DATA 60 


>3. 41 


-0. 9957 


DATA 51 


-3. 33 


-0. 9050 


DATA62 


-1. 42 


-0. 4922 


DATA63 


-2. 41 


-0. 9389 


DAT AS 4 


-1. 03 


-0. 5784 



A simplistic model for sediment volume reverberation was developed and is 
described in Appendix A. Combining the experimental results and the newly developed 
model, the volume scattering coeOicient obtained was 0.0413 ± 0.007 m'^ for fine sand 
and 0.0624 ± 0.007 m for aggregate. These value are reasonable compared to the 
difierence between the decay constant obtained from the experiment. They are 
produced as a best fit using the least square method. The correlation coeflicients, 
approximately 0.88, indicate they are quite a good fit. As a first model, the result is 
promising. However, the discrepancy between the experiment and model might suggest 
the model is still very crude. Also mentioned in the Appendix A is that this is only a 
first step to tr\' to obtain a theory to describe the volume scattering elTect of the 
sediment; it need to be improved by further elTort. 

The ka value for fine sand is 0.12, which is quite small compared to 1. From 
Equation 2.2, cloud-of-small-spheres model, the volume scattering coefficient for the 
fine sand is 0.0082 m''^ and is not too far from 0.0413 m'"’. As mentioned in Chapter 
2, Equation 2.2 was derived under various assumptions and does not apply to real 
sediment. From Equation 2.3, single-small-sphere model, the volume scattering 
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TABLE 6 

AVER,‘\GED nXE SAND REFLECTION 



DATA 


SLOPE 


CORRELATI' 


SET 


( xlO“2) 


COEFFICi: 


DAT 12 


-1. 34 


-0. 9265 


DAT 13 


-2. 09 


-0. 9051 


DAT 14 


-1. 46 


-0. 9309 


DAT 15 


-1. 47 


-0. 9223 


DAT 16 


-2. 36 


-0. 9912 


DAT 17 


-2. 39 


-0. 9904 


DAT 18 


-2. 39 


-0. 9902 


DAT 19 


-2. 38 


-0. 9903 


DAT20 


-2. 37 


-0. 9916 



coefi'ident is 0.34 m''’ which is an order of magnitude different from the the measured 
value 0.0413 m’’^. These values suggest that the cloud-of-small-sphercs model is closer 
to the newly developed model than the single-small-sphere model. 

The scattering coelTicient obtained from Equation 2.2 is 9456 m'^ and from 
Equation 2.3 is 1134 nT'^. Obviously the two values deviate considerably from the true 
value. The ka value is equal to 2.1, which is outside the Rayleigh scattering limit. 

Thus Equations 2.2 and 2.3 can not be applied to the aggregate. 
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VI. CONCLUSION AND RECOMMENDATION 



The following conclusions are possible: 

1. The inhomogenities within the sediment caused the rellected waveform to 
change considerably for various transducer positions. The changes can be 
removed by averaging over the inhomogeneities and the averaging yields a 
mean decay constant. 

2. The glass tank works veiy well, and therefore a bigger tank to contain the 
sediment is not necessary. 

3. The laboratory.' procedure to measure the volume reverberation, using 
off-the-shelf equipm.ent and a high pass filter in a water-over-sediment 
system, is workable. 

4. The basic program, written for a HP-86 desktop computer, is sunicient for 
this experiment. 

5. Although Equations 2.2 and 2.3 can not be applied to aggregate since it 
has a ka greater than 1, the aggregate does have a higher volume 
reverberation than fine sand. The aggregate has a larger size than the 
fine sand; this inight suggest that the larger size sediment causes more 

\ olume reverberation. 

6. From the scattering coefficient 0.0413 ± 0.007 m‘“ for the fine sand and 
0.0624 ± 0.007 m’" for the aggregate, the scattered energy received by 
transducer is quite small compared to the incident energy. 

7. For fine sand sediment, the cloud-of-small-spheres model works better 
than the single-small-sphere model. 

8. The simplistic model, derived from the true acoustical pulse and described 
in Appendix A, gives a crude formulated volume scattering relationship. 
Combining with the experiment result, the sediment's volume scattering 
coefficient can be estimated. 

Further study of the volume reverberation from ocean sediment is recommended. 
Future experimentation should include further variation of the sediment size and of the 
acoustic frequency to find a relationship combining the ka value and volume 
reverberation. The developed model should be improved to obtain a better result 
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coincident with the experimental result. The absorption coelTicient of dilTerent 
sediment with the frequency used should be measured to compare the value obtained 
from Urick's formula described in appendix A. Equations 2.2 and 2.3 should he 
improved to obtain the scattering coelTicient. which can be used to check the above 
model. Continued use of the experimental procedure, equipment and computer 
programs, written for this project on the HP-S6, is highly reconunended. 
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APPENDIX A 

DERIVATION OF SCATTERING FROM THE SEDIMENT 

The generated electrical pulse can be described by the step unit function. II. as 
follows: 



I(t) = Iq ( H(t) - ll(t-T) ) (A.l) 

where l(t) is the intensity at time t, 1^ is the intensity amplitude, t is the pulse length. 

0 t < 0 

M(t) = 

1 t ^ 0 

and 

0 t-T < 0 

H(t-t) = 

1 t-t > 0 

Since the transducer is a resonant system with damping, the square electrical pulse was 
transformed into an acoustic pulse with an exponential rise, an approximately flattened 
top and exponential decay. The shape of the acoustic waveform can be described, 
(Sanders and Coppens, 1986), as 

f(t) = (1 - e'"^) H(t) + (1 - e‘“^ ) ll(t-T) (A. 2) 

where t is the delay time and a is the decay constairt. 

Let a transducer, located in the water, send out a sound pulse to a smoothed 
sediment surface. To assure the surface scattering is negligible and the volume 
scattering is only single scattering, reflected pulse received, by the same transducer, 
should have an intensity (Sanders and Coppens, 1986) 

1 =.(J Tj2 Sy e-Ps f(t-2i'c)d^/(r + ^)2 + R. I(t) ) / r^ (A.3) 
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Where T; is the intensitv transmission cocnicicnt. s„ is the volume scatterin 2 
coclFicient, P is the absorption coelficient in the sediment, t is the alTected sediment 
depth. Rj is the intensity rellection coelTicient and r is the distance between the 
transducer and sediment. From above equation, the intensity contributed by scattering 
is 

s^. e'P^ Iq h:t-2^ c)dc (r + q)- r^ (.A.4) 



and the intensity contributed by the reflection is 
= Ri f(t) r-^ 

Ifr > > q. then 

15 = (1\- s,, r" ;jVr^^t- 2 '^ c;dq 



(A.5) 



{A. 6 ) 



where 



n;t-2^'c) = (1 - c) ) H{t-2^ c) - (1 - ) H(t-t-2^'c) (A. 7) 



During the integration, the lower boundary' of t is zero: the upper boundary- is ct 2 for 
H(t-2i c) and c(t-r) 2 for H(t-r-2i c). Note that only the decay portion, t > t. was 
interested. After the integration, 

I 5 = A ( -(e'^P^'2 - l),p - (e'^P^'2- e‘"V(-P + 2a c) (A. 8 ) 

+ - 1 ) p + );{-p + 2 a/c)) 

where 

A = 1^2 s,. / r^ (A.9) 

From the e.xperimental formula P = 0.25f (Urick, 1979) and fis 182 kHz, the absorption 
coelTicient 
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P = 45.5 dB m = 5.2 Np m. 



cP 2 is 4030 Np s for aggregate and is 41S6 Np s for the fine sand. The intensity decay 
constant a is twice the pressure decay constant 6 . From the experiment. 6 is 37700 
Np s and it gives an a equal to 75400 Np s; 2a c is 97.2 Np m for the aggregate and is 
93.6 Np m for the fine sand. After plugging in the above values for the aggregate. 



Equation .A . 8 becomes 




Ig = A(e“^^^*^^ 11.6) - A(8.3 
For the fine sand. Equation A. 7 becomes 


(A. 10) 


I 5 = A(e“'^^^^ 11.0) - A( 8.6 

R- = R- 
Ri - T^ = 1 


(A. 11) 



The pressure reflection coeflicient R is 0.356 for the aggregate (Diaz, 19S6) and 0.36 
for the fine sand (Bradshaw, 1981). For the aggregate, 

Rj = 0.1270 
Tj = 0.8730 

For the fine sand, 

R- = 0.1296 
T- = 0.8704 

The pulse length is 88 ps and from equation A. 5, 



Ij. = Iq R- (e‘^^ -1) 

Using the above values, the rellected intensity for the aggregate becomes 


(A. 12) 


Ij. = 96.6 Ig /r"^ 

and for the fine sand 


(A. 13) 


= 98.5 Iq /r"^ 


(A. 14) 
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The intensity contributed by the scattering for the aggregate is 



I 3 = (s^. Iq r*^ )( 15.2 - 6.3 ) 



(.A. 15) 



and for the fine sand. 



I 3 = (s^, !q r-^ )( 14.5 - 6.5 ) 



(A. 16) 



Combining the above equations, the reflected intensity from the aggregate received by 
the transducer is: 



I = do B 



(A. 17) 



where 



B = ( s^. 15.2 - 6.3 s^. + 96.6 



400t 



(.A. IS; 



For the fine sand, 

I = (Io'r^)D (A.19) 

where 

D = ( s^. 14.5 - 6.5 s^. + 93 5 g-75400t ^ ^^,^ 20) 

By inputting a series of scattering coefficients and computing the Log (B) *■ and 
1 9 

Log (D) , the theoretical envelope was determined. A least-square-fit slope was 

obtained for each envelope. The scattering coefficient for each sediment can be 

estimated by the value of s„ for which the slope of the theoretical envelope is equal to 

’ » 

the slope of experiment data. For the fine sand, measured the decay constant 23700 
Np;s gave a scattering coefficient of 0.0413 m'^. For the aggregate, the scattering 
coefficient 0.0624 m was obtained from the measured decay constant 21300 Np s. 

Figure A.l and Figure A. 2 show the envelopes determined theoretically from the 
above values of s^, for the aggregate and fine sand respectively. The theoretical 
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envelope is plotted by normalizing the left end of the best fit slope line to the same 
value obtained from the experimient and then applying this normalization to the tvhole 
theoretical envelope. 

Results from the fine sand were selected to find the uncertainty of scattering 
coefficient obtained from* the above method. By comparing a series of e.nveicpe plots 
for different scattering coeificient. the uncertainty = 0.00” was determined. Figure A. 3 
and A.- show the envelopes for scattering coeiTicients w 0.00" diflerent from 0.0^13 
m‘'V The slope changes as the scattering coefficient changes. In the beginning of 
decay, the dominant term in equation .A. 20 is the refection term which is independent 
on s,-. and the scattering term which denends on s,. is neelistibie. The theoretical value 
will not change much in the beginning of decay as the scattering coemcient changes. 
The nonnaiization. described previously, still gave a good approximation to the plot of 
the theoretical envelope vs. the experimental envelope for dhTerent scattering 
coefficients. 

Note that the model derived and described in this appendix is crude. A later 
im^ro cn i.rixS ni.o^ci is 
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LN(-iV0LTAGE'2>'‘l/2) 




Figure A. I For the aggregate, the decay envelope obtained from the model 
vs. the decay envelope obtained fioin the experiment. 
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TIME CMICRC-SED 




Figure A. 2 For the fine sand, the decay enyelope obtained Iroin the model 
vs. tlie decay envelope obtained from the experiment. 



63 



TIME (MICRO-SED 



LN(<V0LTAGE'2>'l/2) 




Figure A. 3 I'or the (ine sand, the decay envelope obtained from the model 
vs. the decay envelope obtained from the e.xpcriment. 
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SLOPE 



TIME CMICRO-SED 



LNC<YDLTAGE*2>*l/2) 
til cL ^ 




Figure A. 4 For the Fine sand, the decay envelope obtained from the model 
vs. the decay envelope obtained IVoin the experiment. 
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APPENDIX D 

COMPUTER PROGRAM 



1 0E>0 
1 0 1 0 
1020 
1030 
1040 
1050 
1060 
1070 
1000 

I 090 
1 1 00 
1110 
1 120 

I I 30 
1 140 
1 150 
1 1 60 
1170 
1 1 90 
1181 
1182 
2000 
20 1 0 
2020 
202 1 
2022 

2023 

2024 

2025 

2026 
2027 
2020 

2029 

2030 

203 1 

2032 

2033 

2034 

2035 

2036 

2037 
2040 
2060 
2061 
2062 

2063 

2064 
2080 
2100 
2120 
2140 

2219 

2220 



I I I I I I ! I i I I I I I I t I I I I I i I I I I t I I I . I I I I I I I I I I t I I ! I I I I I I I I I I I I I > I I I I I I 



PROGRAM : THESIS3 



THIS PROGRAM IS DESIGNED TO CONVERT AND NORMALIZE THE 
THE WAVEFORM DATA VALOE AND POINT LOCATION, TRANSFERED 
FROM THE NICOLET 3091 DIGITAL OSC ILLASCOPE , TO THE 
VOLTAGE AND TIME VALUE. THE NICOLET IS COTROLED BY 
THE HP-3041 DATA ACQUISITION AND CONTROL UNIT. 

THE PROGRAM SQUARE THE VOLTAGE VALUE AND THEN TAKE 
THE RUNNING AVERAGE FOR FIFTY SET OF READING. TO 
INSURE THE DATA IS NOT LOSSED DURING COMPUTER ERROR, 
THE DATA IS BEING STORE IN THE DISK AFTER EVERY 10 
SET CALCULATION. 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I t t I I I t I j I I I I I I I I I I I I I I I • I t ^ I I I 



DIM D(1501) ,T(1501) ,3(1501) 

DIM VDATATC20: , VTIMEFC203 , VC0UNTPC203 
INTEGER N,I 

I 

INPUT THE FILE NAME TO STORE DATA 



DISP "WHAT IS THE FILE NAME FOR VOLTAGE DATA 
INPUT VDATAT 

DISP “WHAT IS THE FILE NAME FOR TIME DATA ?“ 

INPUT VTIMET 

DISP “WHAT IS THE FILE NAME OF THE COUNT DATA ?“ 

INPUT VCOUNT^: 

DISP “YOUR INTEREST TIME INTERVAL IS FROM C START 3 TO CENDT3 “ 
DISP "WHAT IS THE STARTING TIME OF MEASURING IN MICRO-SEC) ?“ 
INPUT START 

DISP “WHAT IS THE ENDING TIME OF MEASURING ( IN MICRO-SEC) ?“ 
INPUT ENDT 



) CHECKING Tf-IE INTERESTED TIME INTERVAL FROM CRT SCREEN 

I 

LET N=1 

I 

IF THIS IS NOT THE FIRST SET OF DATA, THE COMPUTER IS GOING 
! TO READ THE DISK FOR PREVIOUS DAT INFORMATION, AND CARRY ON 
! THE JOB 

I 

DISP “IS THIS YOUR FIRST DATA SET OUT OF THE 50 DATA SET?“ 

DISP “PRESS O IF YES ; PRESS 1 IF N0“ 

INPUT ANS 

IF ANS=1 THEN GOSUB READ^DATA 

I 

! FOLLOWING IS THE COMPUTER RECEIVE THE DATA FROM 3091 AND 
! THEN CALCULATE THE VOLTAGE AND TIME VALUE. 



2240 

2260 

2280 

2290 



DIM ITC53,DTC 20000 3 , NT C 40 3 
CONTROL 10,3 ; 15 

CONTROL 10,4 ; 58 
CONTROL 10,2 ; 1 
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2300 

2301 

2302 

2303 

2304 

2305 

2306 

2307 
230B 

2310 

2311 

2312 

2313 

2314 

2315 

2316 

2320 

2321 



□ UTFUT 709 ;"Q(=N4” 

I 

' ABOVE statement IS LET THE CRT TO "LIVE” CONDITION 

t 

I 

! FOLLOWING ONE IS TO LET THE CRT W "STORE” COrJDITIOM 

I 

I 

OUTPUT 70? ; ”CLS5" 

WAIT 1000 

j 

! FRESS bottom only used without the HP-3421 UTJ I T 
! DISP "FRESS 3091 BUTTON” 

I 

! LET THE DATA TRASMIT TO THE COMPUTER THRDUGHR RS-232 

I 

OUTPUT 709 ;”CLS2” 




ABOVE STATEMENT LET THE WAVEFORM INFORMATION BE TRANSFERED 
FROM CRT TO CCMFUTER THROUGH RS-232 INTERFACE 



2340 ENTER 10 ; It 
2360 ENTER 10 ; Dr 
23B0 ENTER 10 : Nr 
2390 ' 

2400 HO=VAL (NtC 16,203) 

2420 V1^(VAL (NtC21 ,253)-5) (VAL ( N t C 26 . 30 3 ) - 1 2) 

24 4 0 HI = (VAL ( Nt C 3 1 , 35 3 ) -5 ) * 10'" (VAL (Nt C 36 , 40 3 ) - 1 2 ) 

2441 ' 

2443 OUTPUT 709 ; ”0PN5" 

2445 OUTPUT 709 ; ”0PN2” 

2446 ! LET THE CRT IN LIVE MODE 

2447 OUTPUT 709 : ”CLS4” 

2460 J=1 

2490 K=5 

2500 PRINT "N=",N 
2520 FOR 1=0 TO ENDT 

2525 IF P' START OR I , ENDT THEN 2700 

2530 ! 

2531 ! NORMALIZE THE VOLTAGE VALUE AND TIME FOR THE INTEREST INTERVAL 

2532 ! 

2540 D(I)=VAL (DtCJ,K3)*Vl 

2560 T ( I ) =( I-HO) >H1 

25B0 D(I)=D(I) 2 

2600 IF N==l THEN GOTO 2640 

2620 S ( I ) = ( (N- 1 ) *S ( I ) -»-D ( I ) ) /N 

2621 ! 

2640 IF N=1 THEN S ( I ) =D ( I ) 

2660 IF I=START THEM 2760 ELSE 26B0 
26B0 IF I«START+20 THEN 2760 ELSE 2700 
2700 IF I=START*-30 THEN 2760 ELSE 2720 
2720 IF I=START-*-50 THEN 2760 ELSE 2740 
2740 IF I=ENDT THEN 2760 ELSE 2700 
2760 DISP T(I);”SEC ”;D(I>:”VOLT ”;S(I) 

2770 PRINT T(I);”SEC ”:D(I);”VOLT " ; S ( I ) 

2771 ! 

2780 J=J+5 

2B00 M=K. + 5 
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2920 MEXT I 
2040 N = N4-i 
2950 ' 

2051 • LET THE DATA BEING STORE AFTER EVERY 10 SET OF CALCULATION 

2852 ' 

2860 IF N MOD 10=^1 THEM 2?65 ELSE 2240 
2965 IF ^J=ll IHEN 2970 ELSE 29G0 

2970 CREATE VDATAr,200,8 

2971 CREATE VCOUNTr,l ,0 

2972 CREATE VTIMEf,200,8 
2980 GOSUP WRITE_DISK 

2990 IF N=51 THEN 3000 ELSE 2240 
30»:>0 END 

3305 ! 

3306 ! 

3307 ! 

3308 ! FOLLOWING SU&ROUTINE IS USED TO WRITE THE DATA TO THE DISK 

33L)9 ! 

3310 f 

3500 WRITE_DISK: 

3510 PURGE VDATAT' 

3520 CREATE VDA T A.« , 200 , 8 

3540 ASSIGTiH 1 TO VDATAT 

3550 PURGE VCOUMl r 
3560 CREATE VC0UrJTT,l,8 
3580 ASSIGNH 2 TO VC0UT4T r 

3590 PURGE VTIMET 
3600 CREATE V T I ME T , 20 0 , Q 
3620 ASSIGNS 3 TO VTIMET 

3640 PR I NTH 2 ; M 
3660 FOR I=START TO ENDT 
3680 PRIMTtt 3 ; T(I) 

3700 PRINTS 1 : S(I) 

3720 NEXT I 

3740 ASSIGNS 1 TO » 

3760 ASSIGNH 2 TO 
3780 ASSIGNH 3 TO • 

3000 RETURN 

3900 ! 

39 1 0 ! 

3920 ' FOLLOW TNG SUPROUTINE IS USED 10 READ THE DATA FROM PREVIOUS 

3930 ' WRITTEN D I S^ 

394(.) ! 

3950 ' 

4000 READ_DATA: 

4020 ASSIGNH 2 TO VCOUNTT 

4040 ASSIGNH 1 TO VDATAT 

4060 READH 2 ; N 

4061 PRINT N 

4080 FOR I=START TO ENDT 
4100 READH 1 3 S(I) 

4120 IF 1=START THEN 4220 ELSE 4140 
4140 IF I=START+20 THEN 4220 ELSE 4160 

4160 IF I=START+30 THEN 4220 ELSE 4180 

4180 IF I=START+50 THEN 4220 ELSE 4200 

4200 IF I=ENDT THEN 4220 ELSE 4240 

4220 PRINT S(I) 

4240 NEXT I 

4260 ASSIGNH 1 TO 

4280 ASSIGNH 2 TO * 

4380 RETURN 
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801 

802 

8C«3 

804 

805 

806 

807 

808 

809 

810 

811 

812 

813 

814 

815 

316 

880 

900 

910 

'=>20 

930 

940 

950 

960 

970 

980 

990 

991 

992 

993 

994 

996 

1 OOO 

1 020 

1040 

1 060 

1 080 

1 100 

1 120 

1140 

1 141 

1144 

1146 

1148 

1 150 

1152 

1 154 

1 1 60 

1 165 

1 170 

1 180 

1 200 

1220 

1 240 

1260 

1 280 

1290 

1291 



I I I 1 I I [ I I I I I I I I I I I I t 1 I I I I I t I t I I I I I I I I I I « I I I I I I I 1 1 I I I I I I • 1 I I I I 



PROGRAM : MAXIMUM 



THIS PROGRAM IS DESIGNED TO READ THE SQUARED 
WAVEFORM VOLTAGE DATA AND TIME DATA STORED 
IN THE DISK AND THEN CHOOSE THE LOCAL MAXIMUM 
(WAVEFORM ENSEMBLE). THE CHOOSED MAXIMUMS WILL 
BE TAKEN SUARE-ROOT AND STORE THE FINAL RESULT 
ON THE DISK 



I I I I I I I I I I I I I 1 I I « I I I I I I < I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



DIM VDATATC20] ,VTIMETC203 . VMUMBERT C 20 3 , VT I ME _MAXT C 20 3 , VMAX I MUMST C 20 3 
DISP "WHAT IS THE NAME FOR VOLTAGE DATA FILE 
INPUT VDATAT 

DISP "WHAT IS THE NAME FOR TIME DATA FILE 
INPUT VTIMEf 

DISP "WHAT IS THE NAME FOR FILE TO STORE THE NUMBER? NO. ?" 

INPUT VNUMBER? 

DISP "WHAT IS THE NAME FOR FILE TO STORE THE TIME_MAXT DATA ?" 

INPUT VTIME_MAXT 

DISP "WHAT IS THE NAME FOR FILE TO STORE THE MAXIMUMST?" 

INPUT VMAXIMUMST 

DISP "WHAT IS THE STARTING TIME DF YOUR DATA COLLECING ?" 

INPUT START 

DISP "WHAT IS THE ENDING TIME OF THE DATA COLLECTING ?" 

INPUT ENDT 



ASSIGNS* 1 TO VDATAT 
ASSIGN# 3 TO VTIMEr 
CREATE VNUMBERT, 1 ,8 
ASSIGN# 4 TO VNUMBERT 
CREATE VTIME_MAXT, 100,8 
ASSIGN# 5 TO VTIME_MAXT 
CREATE VMAXIMUnST, 100,8 
ASSIGN# 6 TO VMAXIMUMST 



INTEGER FRINTTER 

DISP "WHAT IS THE PRINTER DEVICE ADDRESS MO. ? 
! INPUT THE PRINTER NO. SCREEN IS 1 
INPUT FRINTTER 
PRINTER IS FRINTTER 



DIM 5(1501) ,T(1501) ,TIM(100) ,MAXS(100) 

PRINT USING 1170 ; "TIME ( SEC VOLT " 

IMAGE 25X , 1 lA, 12X ,4A 
INTEGER I,J 
FOR I==START TO ENDT 
READ# 3 ; T(I) 

READ# 1 ; S(I> 

NEXT I 
LET J=^0 

i 

! IF THE FIRST DATA BEING READ IS GREAT THAN THE SECOND ONE, 
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1292 

1I?3 

1 300 

1320 

1340 

1360 

1380 

1 40<:j 

1410 

1420 

1440 

1460 

1480 

1500 

1520 

1540 

1560 

1580 

1 600 

1610 

1620 

1640 

1660 

1680 

1 700 

1720 

1730 

1740 

1760 

1 800 

1820 

1840 

1 880 

1 900 



! THE FIRST ONE IS TREATED AS A LOCAL riAXIMUN. 
1 

IF S(START>:-= S(START>1) THEN 1320 ELSE 1500 
TIM(J)=T(START> 

MAXS(J)=SaR (S( START) ) 

PRINTS 5 j TIM(J) 

PRINT# 6 ; NAXS(J> 

PRINT USING 1410 ; " J= ” ; J , T I M < J ) , MAXS ( J ) 

IMAGE 12X ,3A,DDD,7X,D. DDDDDD, 1 IX ,D. DDDDDDDD 
! PRINT TIM(J) 

! PRINT MAXS(J) 



FOR I=START+1 TO ENDT-1 

IF S(I>. = S(I-l) AND S(I).;= S(I-*"1) THEN 1540 ELSE 1720 
TIM(J)=T(I) 

MAXS ( J) =SDR (S ( I ) ) 

I 

PRINT USING 1610 ; ” , J , T I M ( J ) , MAXS ( J ) 

IMAGE 12X ,3A,DDD,7X ,D. DDDDDD, 1 1 X , D . DDDDDDDD 
! PRINT TIM(J) 

' PRINT MAX3(J> 



PRINT# 


cr - 

^ ! 


TIM(J) 


PRINT# 


6 ; 


MAXS ( J) 


J=J + 1 








NEXT I 








J=J-1 








PRINT# 


4 ; 


J 




ASSIGN# 


4 


TO 




ASSIGN# 


5 


TO 


♦ 


ASSIGN# 


6 


TO 




ASSIGN# 


1 


TO 


M- 


ASSIGN# 


3 


TO 




END 
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90^1' ' ' I > j I I I I I I I I I I I I I I I I I I I ) 1 I I I I I I I I I I I ' I I I I I I I I 1 I I I I I I I I I I I I < I I 

^01 ! ' ! 

903 ! ! ' 

'=’04 ! ! PROGRAM ; PL0T3 1 

905 ' ! ! 

906 ! ' ! 

907 ! ! THIS PROGRAM IS DESIGNED TO PLOT THE ENVELOP ! 

90S ! ! OP WAVEFORM BY READING DATA FROM DATA DISK. ! 

909 ! ! AFTER READING, VOLTAGE DATA IS CONVERTED TO LOG ! 

"=^10 ! ! VALUE, THEM BE PLOTTED. AFTER PLOTING, A LEAST- ! 

911 ! ! SQUARE METHOD IS USED TO FIND THE BEST FIT STRAIGHT ! 

912 ! ' LINE FOR THE LOG-ENSEMBLE, PLOT THE LINE, PRINT ' 

913 ! ! OUT THE SLOPE OF THE LITJE AND THE CORRELATION ' 

914 ! ! COEFFICIENT- , ! 

915 ! ! ! 

916 • ' ! 

917 I I I I I I I I I I I I I I I I • I I I I I I I I I I I I I 1 I I I I I I I I I 1 I I I I I 1 1 I I I I I I I I I 1 I I 1 I 

918 ! 

919 ' 

1000 CLEAR 
1020 GCLEAP 

1022 INTEGER PLOTER 

1024 DISP ’‘WHAT IS YOUR POLTTER DE’/'ICE ADDRESS NO, ?'* 

1025 ! INPUT THE PLOTTE NO. SCREEN IS 1 

1026 INPUT PLOTER 
1028 PLOTTER IS PLOTER 

1040 FRAME 

1041 ! 

1042 ! THIS PLOT SIZE IS TO FIT THE REQUIREMENT FOR THE NFS THESIS 

1043 ! 

1060 LOCATE 30,125,30,85 

1080 CSIZE 3 

1100 MOVE 50,20 

1120 LABEL “TIME (MICRO-SEC)" 

1140 PEN UP 
1160 MOVE 18,40 
1180 DEG 
1200 LDIR 90 

1220 LABEL "LN ((: VOLTAGE 2> ) 1 /2 ) " 

1240 PEN UP 

1241 ! 

1242 ! DEFAULT SET-UP IS TO FIX THE PLOTING COORDINATION AS 

1243 ! X-AXIS IS FROM 500 TO 670 (MICRO-SEC), Y-AXIS IS FROM 

1244 ! -5 TO 2 

1245 ! 

1250 DISP “DO YOU WANT TO USE DEFAULT SET-UP YES-0 “ 

1251 INPUT BULL 

1252 IF BULL=0 THEN 1253 ELSE 1260 

1253 SCALE 500,670,-5,2 

1255 LAXES 1 0 , . 5 , 500 , — 5 , 2 , 2 

1256 XMIN-500 

1257 GOTO 1710 

1259 ! 

1260 DISP “ENTER THE XMIN OF SCALE" 

1280 INPUT XMIN 

1300 DISP “ENTER XMAX OF SCALE" 

1320 INPUT XMAX 

1340 DISP “ENTER THE YMIN OF SCALE" 

1360 INPUT YMIN 
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1380 

1 400 

1420 

1440 

1460 

1480 

1 500 

1520 

1540 

1560 

1 580 

1600 

1620 

1640 

1 660 

1661 

1 680 

1700 

1710 

1720 

1721 

1740 

1741 

1742 

1744 

1746 

1748 

1749 

1750 

1751 

1752 

1759 

1760 

1780 

1781 

1782 

1790 

1800 

1802 

1803 

1 804 

1810 

1820 

1821 

1822 

1823 

1825 

1826 

1828 

1831 

1832 

1833 

1834 

1 835 

1837 

1839 

1840 

1 920 

1930 

1940 



DIBP "ENTER THE YMAX OF SCALE 
INPUT YMAX 

SCALE XMIN,XNAX,YMIN,YMAX 
FXD 0,2 



DI3F "ENTER 
INPUT XT 
DISP "ENTER 
INPUT YT 
DISP "ENTER 
INPUT XI 
DISP "ENTER 
INPUT YI 
DISP "ENTER 
INPUT XMC 
DISP "ENTER 



THE X-TICIING SPACE" 
THE Y-TICriNG SPACE" 
THE X INTERSECTION" 
THE / INTER SECTION" 
THE X-MAJOR COUNT" 
THE Y-MAJOR COUNT" 



INPUT YMC 

LAXES XT,YT,XI ,YI ,XNC,YMC 
GRAPH 
FEN 1 



DIM TIM ( 100) ,MAXS ( 100) 

Din TIMEMAXTCCO] . MA X I MUMT C 20 ] 

DISP "WHAT IS YOUR FILE NAME FOP TIMEMAXT 
INPUT TIMEMAXT 

DISP "WHAT IS YOUR FILE NAME FOR MAXIMUM^ ?" 

INPUT MAXIMUM^ 

PRINT "THE DATA FILES APE " , T I MEMAXT , " AND " , MAX I MUMT 
DISP "WHAT IS THE HIGHEST ORDER OF ARRAY NUMBER';'" 

! THIS GIVES A CHODICE TO CHOOSE THE END POINT OF PLOT 
INPUT JK 

I 

ASSIGNl^ 5 TO TIMEMAXT 
ASSIGN# 6 TO MAXIMUM^ 

DISP "PRESS O IF THE FIRST LOCAL MAXIMUM IS BEING SAVED ?" 

INPUT OK 

FOR 1=0 TO JK 

READ# 5 ; TIM(I) 

! CONVERT THE TIME UNIT FROM SEC TO MICRO^SEC 

TIM(I)=TIM(I)^1 000000 
READ# 6 ; MAXS(I) 

MAXS(I)=LOG (MAXS(I) ) 

! CONVERT THE VOLTAGE TO NATURE LOG VALUE 

NEXT I 

K=JK+1 

IF 0K=0 THEN GOTO 1835 
FOR 1=1 TO JK 
TIM ( I“1 ) =TIM ( I ) 

MAXS(I-1)=MAXS(I) 

NEXT I 
K=JK 

! DISP "WHICH CHARACTER YOU PREFER TO USING IN PLOTTING THE DATA? 

! INPUT CT 

MOVE TIM (0) ,MAXS (0) 

FOR 1=0 TO K-1 
DRAW TIM ( I ) ,MAXS ( I ) 

! LABEL CT 
NEXT I 
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1960 

1961 

1962 
19G0 
2000 
2020 
2900 
2910 
2920 
2922 
2924 
2926 
2928 
2930 
3000 
3040 
3060 
3080 
3100 
3120 
3140 
3160 
3 1 SO 
3200 
3220 
3240 
3260 
3280 
3300 
3320 
3340 

3350 

335 1 



3360 

3370 

337 1 

3372 
3380 
3400 

3420 

3421 

3422 

3423 

3440 

3441 

3442 

3443 

3460 

3461 

3480 

3481 
3500 

3520 

3521 
3540 

3560 

3561 

3562 
3580 



FEN UP 

GOSUB LEASTSQR 
ASSIGN^# 5 TO 
ASSIGN# d TO •* 
END 



‘ FOLLOWING SUBROUTINE IS DESIGNED TO APPLY THE LEAST 

! SQUARE METHOD IN THE ABOVE DATA AND THE PLOT THE 

! BEST FIT STRAIGHT LINE, PRINT OUT THE SLOPE AND 

! THE CORRELATION COEFFICIENT ON THE PLOT. 

! FORMULA USED PLEASE 

1 

LEASTSQR: 

REAL Y , X Y , X 1 , X2 , TAVE , YAVE , VARX , VARY , COVX Y 
LET Y=0 

! Y IS THE SUMMATION OF Yd) VALUE 

LET XY=0 

! XY IS THE SUMMATION OF X (I ) *Y ( I) 

LET X1=0 

! XI IS THE SUMMATION OF X(I) 

LET X2=0 

! X2 IS THE SUMMATION OF X(I) 2 

FOR 1=0 TO K-1 
Y=Y+MAXS ( I ) 

XY=XY+TIM ( I ) ♦MAXS (I) 

X1=X1+TIM(I) 

X2=X2+TIM(I) 2 
NEXT I 

B= (X2^Y-X 1*XY) / (K>X2-X1 '2) 



! A IS THE SLOPE OF THE FITTED STRAIGHT LIGHT 

A=(P>XY-X1*Y) / (K*X2-X1 2) 



! B IS THE INTERCEPT TO Y AXIS 

PRINT "THE SLOPE A=";A 
PRINT "THE Y- INTERCEPT B=";B 
YAVE= (MAXS (0) +MAXS (K-1 ) ) /2 

I 

! YAVE IS THE AVERAGE VALUE FOR THE LOG VOLTAGE 

I 

TAVE=(TIM(0)+TIM(K-1) ) /2 

I 

? TAVE IS THE AVERAGE VALUE FOR THE TIME VARIABLE 

I 

C0VXY=0 

i 

! COVXY IS THE COVARIANCE OF TIME AND VOLTAGE (LN) 

I 

VARX=0 

! VARX IS THE VARIANCE OF TIME VARIABLE 

( 

VARY=0 

! VARY IS THE VARIANCE OF VOLTAGE (LN) VARIABLE 

DIM YL(IOO) 

FOR 1=0 TO K-1 
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3600 

3620 

3640 

3 6 6 

36S0 

3700 

3720 

3740 

3760 

3761 

3730 

3800 

3820 

3840 

3860 

3880 

3900 

3910 

3920 

3940 

3960 

3980 

3990 

3991 

3994 

3996 

3998 

3999 

4000 

4002 

4004 

4006 

4020 

4021 

4040 

4060 

4061 

4062 

4080 



C0VXY=C0VXY'»' (TIIi< I ) -TAVE) ^ (MAXS ( I ) -YAVE) 
» 

VARX = VARX+ (TIM ( I ) -TAVE) ^-^2 

i 

VARY=VARY+ (MAXS ( I ) -YAVE) 



YL(I)=A*TIM(I)+B 

! YL IS THE VALUE FITTED ON THE STRAIGHT LINE 



NEXT I 
PEN 2 

MOVE TIM (0) , YL (O) 

FOR 1=0 TO K>1 
DRAW TIM(I) ,YL(I) 

NEXT I 

CORCOE = COVXY/SDR (VARX-^VARY) 

I 

! CORCOE IS THE CORRELATION COEFFICIENT 

I 

PRINT "CORRELATION COEFFICIENT IS CORCOE 

PEN I 

DEG 

LDIR 0 

I 

! FOLLOWING PRINT OUT THE SLOPE AND COORELATION COEFFICIENT 
! ON THE PLOT 

I 

MOVE XMIN+80,1 

I 

! LABLE ALWAYS STARTS FROM THE XMIN-^80 IN X COORDINATES 

LABEL USING 4021 ; "SLOPE =",A 
IMAGE 7A,DDD.DDDD 
MOVE XMIN+80,2 

LABEL USING 4061 ; "CORRELATI OIJ COEFF= ", CORCOE 

IMAGE 19A,DDD. DDDD 
PEN UP 
RETURN 
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BOO 

010 

S20 

sro 

S40 

850 

860 

070 

880 

090 



I 1 ) I I I I I I I I I I I I I I i I I i I I I 1 I I I t I I t I I I I I I I I I I I t I I < I 



PROGRAM ; AVERAGE ! 

THIS PROGRAM IS DESIGNED TO AVERAGE ! 

THE DATA SETS TO GET A AVERAGED ! 

WAVEFORM DATA. ! 

I I I I • I I I I I I I I I I I I I I I I I I I I I I I I I ( I I I I I I I I I I I I I I I I I I 



1000 DIM VDAT^C253 , VDATA^ C253 , VDATU^C253 
1020 DIM VI ( 1500) ,V2( 1500) ,V3( 1500) 

1040 DISP "WHAT IS THE ORIGINAL DATA FILE ? DAT#" 

1060 INPUT VDATT: 

1080 DISP "WHAT IS THE NEW DATA FILE ? DATA#" 

1100 INPUT VDATAT^ 

1120 DISP "WHAT IS THE AVERAGED DATA FILE ? DATl#" 

1140 INPUT VDATIT: 

1160 ASSIGN# 1 TO VDATJ^ 

1180 ASSIGN# 2 TO VDATAT 

1200 CREATE VDAT 1 4= , 8 , 1 10 
1220 ASSIGN# 3 TO VDATIT 

1 230 ! 

1231 ! N IS THE NUMBER HOW MANY DATA SET BEING AVERGED 

1232 ! 

1240 DISP "WHAT IS YOUR AVERAGE No. N ?" 

1260 INPUT N 

1261 ! 

1262 ! AVERAGE ONLY APPLY TO THE SAME ARRAY NO. 

1263 

1270 FOR 1=522 TO 621 

1271 READ# 1 ; VI (I) 

1272 NEXT I 

1274 FOR 1=522 TO 621 

1275 READ# 2 ; V2 (I ) 

1276 NEXT I 

1280 FOR 1=522 TO 621 

1 300 ! 

1301 ! RUNNING AVERAGE IS ADOPTED 

1302 ! 

1340 V3 ( I ) = ( (N-1 ) ^V1 ( I ) +V2 ( I ) ) /N 
1360 PRINT# 3 ; V3 ( I > 

1371 ! 

1372 f PRINT OUT FOLLOWING VALUE AS A CHECKING PURPOSE 

1373 ! 

1380 IF 1=522 THEN 1500 ELSE 1400 

1400 IF 1=553 THEN 1500 ELSE 1420 

1420 IF 1=600 THEN 1500 ELSE 1440 

1440 IF 1=621 THEN 1500 ELSE 1520 

1500 PRINT VI (I) ,V2(I) ,V3(I) 

1520 NEXT I 

1540 END 
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APPENDIX C 

DECAY ENVELOPE GRAPHS FOR AGGREGATE 

This appendix includes the graphs of the decay envelopes for reflection from 
aggregate water interface within the glass tank. Figure C.l to C.9 show the original 
decay envelopes obtained from different transducer positions. The decay envelopes are 
also described as following data sets: 

data set 42 
data set 43 
data set 44 
data set 45 
data set 46 
data set 47 
data set 
data set 49 
data set 50 

Figures C.IO to C.l 8 arc the decay envelopes resulted from the running average 
method. They are also described as following data sets: 

data set DAT2 
data set DATS 
data set DAT4 
data set DAT5 
data set DAT6 
data set DAT7 
data set DAT8 
data set DAT9 
data set DAT 10 
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LN((<VDLTAGE‘2>) ‘1/2) 




Figure C. 1 Decay envelope for rdlection fVoin aggregate/water interface, 
transducer was 30 cm above surface (data set 42). 
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l igure C.2 Decay envelope for reflection from aceregatc, 'water interface, 
transducer was 3U cm above surface (data set 43). 
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TIME (MICRO-SEC) 
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Figure C.3 



Decay envelope for rellectioii from aggregate, 'water interface, 
transducer was 30 cm above siuiace (data set 4d). 
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CORRELATION C0EFF= -.7848 



TIME (MICRO-SEC) 



LN((<VDLTAGE*I?>) *1/2) 




Figure C.4 



Dccav envelope lor redection from agcrcgatc^vat 
transducer was 30 cm above surface (data set 45). 
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Figure C.5 Decav envelope for reflection from ag^regate/water interface, 
transducer was 30 cm above surface (data set 46). 
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Figure C.6 Decay envelope for relicction from agmegate/water interface, 
transducer was 30 cm above surface (data set 47). 
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CORRELATION C0EFF« -.3971 



LN((<V0LTAGE*2>) *1/2) 




Figure C.7 Decav envelope for rcllcction from aggvegate/watcr interface, 
transtfucer was 30 cm above surface (data set 48). 
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Figure C.8 Dccav envelope for redection from aggregate/water interface, 
transducer was 30 cm above surl'ace (data set 49). 
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CORRELATION C0EFF“ -.8654 




Figure C.9 Dccav envelope for rellection from aggregate/water interface, 
transducer was 30 cni above surface (data set 50). 
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I'igurc C. 10 Decay envelope resulted from the average of data set 42 and 43 
for the aggregate/water surface (data set DAT 2). 
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. OOr CORRELATION C0EFF= -.6051 
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F'igure C.l 1 Dccav envelope resulted (Voin llie average from data setd2 to 44 
lor the aggregate/water surl'ace (data set DA 13). 
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CORRELATION C0EFF= -.9720 



LN((<V0LTAGE'2>) ‘l/'2) 




Figure C.12 Occav envelope resulted (roin the average from data set42 to 45 
for the aggregate/water surface (data set DA'14). 
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Figure C. 13 Dccav envelope resulted from the average from data set 42 to 46 
for the aggregate/ water surl'acc (data set DA 1'5). 
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Figure C.14 



Decay envelope resulted from the average from data set 42 to 47 
for the aggregate, 'water surl'ace (data set DA'Kj). 
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I'igurc C. 15 Dccav envelope resulted from tlie average from data set 42 to 48 
lor the aggregate/water surl'ace (data set DA 17). 
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Figure C.16 Dccav envelope resulted IVoni the average from data set 42 to 49 
(or the aggregate/water surlhce (data set DATS). 



92 



CORRELATION C0EFF“ -.9932 
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Figure C.17 Dccav envelope resulted from the average from data set 42 to 50 
for the aggregate/water surface (data sel DA 19). 
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Figure C.18 Decay envelope resulted from the average from data set 42 to 51 
for the aggregatc/water surface (data set DATIO). 
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APPENDIX D 

DECAY ENVELOPE GRAPHS FOR FINE SAND 



This appendix includes the graphs of the decay envelopes for reflection from fine- 
sand, water interface within the glass tank. Figure D.l to D.IO shew the original decay 
envelopes obtained from different transducer positions. The decay envelopes are also 
described as following data sets; 

data set 55 
data set 56 
data set 57 
data set 58 
data set 59 
data set 60 
data set 61 
data set 62 
data set 63 
data set 64 

Figures D.ll to D.l 9 are the decay envelopes resulted from the running average 
method. They are also described as following data sets: 

data set DAT 12 
data set DAT 13 
data set DAT 14 
data set DAT 15 
data set DAT 16 
data set DAT 17 
data set D.'\T1S 
data set DAT 19 
data set DAT20 
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Figure D.l Decay envelope for reflection from fine-sand/water interface, 
transducer was 30 cm above surface (data set 55). 
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Figure D.2 Decav envelope for rellcction from fine-sand/water interface, 
transducer was 30 cm above surface (data set 56). 
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Figure D.3 Dccav envelope for reflection from fme-sand/water interface, 
transducer was 30 cm above sui face {data set 57). 
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Figure D.4 Decay envelope for rcllcction from nne-sand, water interface, 
transducer was 30 cm above surl'ace (data set 58). 
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Figurc D.5 Dccaj' envelope for rcllcction from fine-sand/water interface, 
transducer was 30 cm above surface (data set 59). 
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Figure D.6 Dccav envelope for reflection from (Ine-sand/water interface, 
transducer was 30 cm above surface (data set 60). 
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Figure D.7 Decay envelope for rdkction from line-sand, 'water interlace, 
transducer was 30 cm above surface (data set 61). 
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Figure D.8 Dccav envelope for rellection from nne-sand water interface, 
transcilicer was 30 cm above surface (data set 62). 
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Figure D.9 Dccav envelope for renection IVom fine-sand, water interface, 
transducer was 30 cm above surface (data set 63). 
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Figure D.IO Decay envelope (or redection I'roin fine-sami/water interface, 
transducer was 30 cm above surface (data set 64). 
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Figure D.l 1 Decay envelope resulted from the average of data set 55 and 56 
I’or the fme-sand/water surl'acc (data set DA'1'12). 
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Figure D.12 Decay envelope resulted from (he average from data set 55 to 57 
for the fme-sand/water surface (data set DA 1 13). 
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Figure D.13 Decav envelope resulted from the average from data set 55 to 58 
for tlie fme-sand/water surface (data set DA'f 14). 
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Figure D.14 Decav envelope resulted from the average from data set 55 to 59 
lor the line-sand/water surface (data set DAT 15). 
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ORRELATION CCEFF= -.9912 
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Figure D. 16 Dccav envelope resulted (Voin tire average from data set 55 to 61 
lor tlie line-sand/water surface (data set DA 117). 
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Figure D.17 Decay envelope resulted froin the average from data set 55 to 62 
lor the line-sand/water surl'ace (data set DAT 18). 
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Figure ID. 18 Dccav envelope resulted lioiu the average Iroiu data set 155 to 63 
Cor the linc-sand/watcr surracc (data set DA 119). 
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I igurc D.19 Dcc.'iv cnvcloj)c icsultcd IVoin llic ;ivcr:igc (Voni data set 55 to 64 
lor tlic aggregate, 'water suilace (data set D/\ 1 20). 
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